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ABSTRACT 

The 0{am^ /m%/) Yukawa corrections to top pair production in photon-photon collision 
are calculated in the standard model (SM), the general two-Higgs-doublet model (2HDM) 
as well as the minimal supersymmetric model (MSSM). We found that the correction to the 
cross section can only reach a few percent in the SM, but can be quite significant (>10%) 
in the 2HDM and MSSM for favorable parameter values, which may be observable at the 
high energy e^e~ colliders. 
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1. Introduction 

Recently, the evidence for top quark production, with a mass of 176 it 8{stat) it 10{syst) 
GeV and 199'^l^{stat) ± 22{syst) GeV has been reported by the CDF and DO collaboration, 
respectively [1]. Due to its large mass, the discovery of the top quark will open a number 
of new and interesting issues, such as the precision measurement of the mass, width and 
Yukawa couplings of the top quark through its direct production and subsequent decay at 
both hadron and e^e~ colliders. But even with 1000p6~^ of luminosity, the Fermilab Tevatron 
could determine the top mass to 5 GeV or better [2]. At the future multi-TeV proton colliders 
such as the CERN Large Hadron Collider (LHC), tt production will be enormously larger 
than the Tevatron rates, but the accuracy with which the top mass can be measured in proton 
colliders is limited to about 2 ~ 3 GeV[2]. Bloude et al.[3], have argued that one must know 
the top mass to 1 GeV to take full advantage of the constraints that precision electroweak 
measurements put on the Higgs boson and other massive particles which might contribute to 
electroweak loops. Beyond this, it would be wonderful to make a precision measurement fo 
the basic parameter nit to 0.3GeV or better for looking for new physics beyond the SM by 
the loop processes which are sensitive to mj. At the next-generation linear collider (NLC) 
operating at a center-of-mass energy of 500 GeV-2000 GeV with a luminosity of the order 
of 10^^ cm^^s"^, the e'^e" — tt events rate would be around 10^/yr, comparable with the 
Tevatron, however the events would be much cleaner and top parameters would be easier 
to extract. At the NLC a top mass measurement with statistical uncertainty 0.3 GeV from 
10/6~^ luminosity is expected [2] and it is possible to separately measure all of the various 
production and decay from factors of the top quark at the level of a few percent [4]. 

Nowadays, the possibility of transforming a linear e+e" collider into a 77 collider deserves 
a lot of attention. With the advent of the new collider technique [5] the collision of high 
energy, high intensity photon beams , obtained by using the old idea of Compton laser 
backscattering[6], can be realized in the NLC. The back-scattering of laser photons off the 
colliding electron and positron beams can yield intense and energetic photon beams which 
then collide with the high luminosity. There are many uses for such high photon-photon 
luminosity, one of the most important may be for the production of top quark paris. It 
has been found [7] that tt production in 77 collisions realized by laser back-scattering is 
slightly larger than the direct e'^e~ — > tt production for mt < 130 GeV at ^/s = O.BTeV 
, and at s/s = ITeV the production of 77 — tt is much larger than e+e~ — tt for mt ~ 
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100 — 200 GeV both with and without considering the threshold QCD effect. In the SM, 
the cross section for top quark pair productions in 77 coUisions have been calculated with 
higher order QCD correction [8]. The radiative corrections to 77^ from final state Higgs 
exchange interactions has also computed in Ref.[9]. The correction is of order 0(2 — 4%) for 
typical values of the Higgs boson mass and top quark mass. In this paper we calculate the 
0{amf/m^r) Yukawa correction in a two Higgs doublet model (2HDM) (Model II)[10] and 
in the minimal supersymmetric model(MSSM) , in which there are three neutral and two 
charged physical Higgs bosons, H, h, A, , of which H and h are CP-even and A is CP- 
odd. The 0{am^ /m'^r) Yukawa correction arise from the virtual effects of the third family 
(top and bottom) quarks, charged and neutral Higgs bosons, as well as the Goldstone bosons 
(G°,G^). The results of the standard model can be obtained from our calculations as a 
special case. In Sec. II, we present the analytic results in terms of the well-known standard 
notation of one-loop Feynman integrals. In Sec. HI, we present some numerical examples 
and discuss the implication of our results. And in the appdendix we list the form factors in 
the cross section. 

2. Calculations 

The relavent Feynman diagrams are shown in Fig.l and the Feynman rules can be found 
in Ref.[10]. In our calculation, we use dimensional regularization to regulate all the ultraviolet 
divergences in the virtual loop corrections and we adopt the on-mass-shell renormalization 
scheme [11]. In our calculations we keep the term mfetan/? in the in the charged Higgs 
couplings to third family quarks since its effects become rather important for large tan/?. 
Taking into account the 0{am^ /m^r) Yukawa corrections, the renormalized amplitude for 
77 tt is given by 

Mren = Uj:^^ + Mj^^l (1) 

= M^*^ + (^M*(*) + (^M''^*) + 5M''(*) + (^M^(*) , (2) 
= M(^^{ps^p,,i^u), (3) 

where Mq is the amplitude at tree level, SM'^ , 5M'" , 5M^ and SM^'^*'^ represent the 0(am^/m^) 
Yukawa corrections arising from the self energy diagram Fig. 1(d), vertex diagram Fig.l(f)-(i) , 
box diagrams Fig.l(l)-(n) and digrams Fig.l(j),(k), respectively, i = (^4—^2)^, u = {Pi—Pa}^ 
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and ^3(^4) denote the momentum of the two incoming photons, and P2(pi) are momentum 
of the outgoing t quark and its antiparticle. The exphcit forms of these matrix elements are 
given by 



e2 



t-rrit 



SM^(') = ie^Ql , e^{p^Y{P3)u{v2){jl 7^7. + !l Pi^ 1, + /e P2m 7. 

+/r2 ;^47m7i/ + /la i>^P\vlp. + /le l^4P2/.7!/)^'(Pi), (6) 

+/Ip5'7'' + fVxi" + /l^>27'^ + flp^i"" + /M^r 

+/IKP2 + /gP^rf + /iV^P2 + /nS'^^l^) + /i2;^47^^ 

+/l'3;^4Pr7'^ + /W^" + /l5;^4P2^7'^ + /l6^^4P^7" 

+/f7/^4P^^rf + A^hP^xPl + /l9«Pl + /lo«i'2] ^(Pl) (7) 

<5M^(*) = ze^Q L 2 g^(P4)e.(P3)n(p2) [/iV^ + /tVK 

+/8Vf'2 + /g^'i'^K + /i''oP^i'2l ^(i^i) (8) 



The form factors /f are presented in Appendix A. 

The corresponding amphtude squared can be written as 

+(5MWM(S'')^) (9) 
is given by 

Y^m^'^'^Mf = 5M'^')Mf + 5M<'^Mf + SM^^'^uf + SM^'^'^Mf , (10) 

where 



^JM^WmW^ = ,, (/|g2 + /|iJ6 + /f2gl2), (11) 

v3^2r)2 



^ 4:m^s^{t - mi) 
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+f^2Hi2 + /rs^is + fwHie), (12) 

^JM^'WMf)' = ^^!pl^Y^f^H,{mt,p,-p2,Pi-P3,Pi-p,,Ps-p,,), (13) 

4m^s^(t-m^)^ 

'im^sf^it - mf) 

Here the expressions of the Hi{mt, pi • P2, Pi • Ps, Pi • Pa, Ps • Pa) are given in Appendix B. 

can be obtained by 



JM^Mo^")^ = (5MWMo(*)t|p^^p^ (15) 
(^M^Mo^^ = (5mWMoW^|p4^p3 (16) 
(5M(")Mo(")^ = (5mWMo(*)^|p4„p3 (17) 

The cross section of the subprocess is given by 

(t(s) = ^^^* di\Mrenis,i)\\ (18) 

where = {mf — ^s) ± ^s^t and (3t = y^l — Amf / s . The total cross section for top quark 
pair production can be obtained by folding the cross section a for the subprocesses with the 
photon luminosity 

rXmax _ 

^{^) — / dzdL^-y/dza^jj tt sit s = z s) (19) 

where ■\/s{-\/§) is the e'*"e~(77) center-of-mass energy and the quantity ^^^^ is the photon 
luminosity, which is defined as 

^=2-/Tf^7/e(-)^./e(-V-) (20) 

^max 

For unpolarized initial electrons and laser, the energy spectrum of the back-scattered photon 
is given by [12] 



where 



1 4x 4x^ 
1-X + - — - + 



1-x c(i-x) ei^-xy 



(21) 



m) = (1 - 7 - 72 ) 1^(1 + + o + 7 - TJTT^ (22) 



1 8 1 
e e'"''' ' ^' ' 2 + 2(1 + 0^ 
and ^ = 4^oWo / ' and £'o are the incident electron mass and energy, respectively, and 
ujQ is the laser- photon energy, x is the fraction of the energy of the incident electron carried 



by the back-scattered photon. In our calculation we follow the analysis of Ref. [9], and choose 
luq such that it maximizes the back-scattered photon energy without spoiling the luminosity 
through e+e~ pair creation. With this choice, we can finds ^ = 2{l + \/2) ~ 4.8, Xmax — 0.83, 
and D(0 ^ 1.8 . 

3. Numerical results and conclusion 

In our numerical calculation, the input parameters [13] arc mz = 91.176GeV, aem = 
1/128.8, and Gf = 1.166372(10-5(GeV)-2). mw is determined through [4] 

where, to order 0{am1 /my^), Ar is given by [4] 

The lower limit of the parameter tan/3 is 0.6 from perturbative bounds [14]. Reference [15] 
argues lower values of tan/3 from perturbative unitarity which is about 0.25 for top quark 
mass of 176 GeV. So in our numerical calculations we allow tan/3 to take the minimum value 
of 0.25 in the two Higgs doublet model. In the following we present some numerical examples 
corresponding to a e^e~ collider with center-of-mass energy of -y/s = 500 GeV. 

The numerical results in the SM are presented in Fig.2. The correction to the cross section 
depends on the Higgs mass and at Mh = 300 GeV the correction reaches its maximum size 
of —2.7%. Recently, the correction in the standard model has been calculated in Ref. [16]. 
But in that work the authors only present the correction to subproccss cross section 77 — > tt 
and did not give the corresponding results at e+e" collider. So it is difficult to compare their 
results with ours. 

We present the numerical results in the two-Higgs-doublet model in Fig.3 and Fig.4. In 
our results we fix the parameters a and /3 to be a = /3 = 0.25 and show the dependence 
on the masses of Higgs bosons. The correction is sensitive to the Higgs masses and can be 
quite large for small Higgs masses. Fig.3 shows the dependence on the mass of CP-even 
Higgs bosons h and H for fixed Ma value. Wc found that correction can be quite large for 
small Mh value. For M^ < 100 GeV the correction can exceed 50% and makes it necessary 
to calculate higher order corrections beyond one-loop level. Fig.4 shows the dependence on 
the mass of CP-odd Higgs boson A for fixed M^^h value. For Ma = 100 GeV the correction 



reaches -38% and decreases rapidly with the increase of Ma_. The corrections drops rapidly 
with the increase of tan /3 as the case of the minimal supersymmetric model discussed bellow. 
Here we did not present the numerical results corresponding to large tan f3. 

Figs.5-7 represent some numerical results in the minimal supersymmetric model. The 
Higgs sector of the minimal supersymmetric model is a special case of the two-Higgs-doublet 
model. In this model the masses and couplings of the Higgs bosons are controlled by two 
parameters at tree level, which can be taken to be Ma and P for example. In our numerical 
results presented in Figs. 3-5, we show the dependence on Ma for three different values of 
tan p. From these figures one can find that the correction depends strongly on the values 
of tan/3. The correction is more significant for smaller tan/3 values. And for a fixed tan/? 
value the correction can be either positive or negative, depending on the Higgs mass Ma- For 
minimum (3 value (3 = 0.25, the correction gets its positive maximum size of 13% at Ma = 420 
GeV and negative maximum size of -54% at Ma = 300 GeV. For tan (3 = 1 the positive and 
negative maximum size of the correction can only reach 7% and -1.6%, respectively. For 
larger tan/3 value tan/3 = 5, the behaviour of plot in Fig. 7 is different from small tan (3 plots 
in Figs. 5, 6 since the effect of the coupling ~ tan (3 becomes significant when tan (3 gets 
large and cancel to some extent the effect of the coupling mt cot f3. 

In conclusion, we calculate the 0{arn^ /m'^Y) Yukawa corrections to top pair production 
in photon-photon collision in the standard model (SM) , the two-Higgs-doublet model as well 
as the minimal supersymmetric model. We found that the correction to the cross section 
can only reach a few percent in the SM, but can be quite significant (>10%) in the 2HDM 
and MSSM for favorable parameter values. So these corrections are potentially observable 
at next generation linear collider, and thus could be used to set limits on the parameters of 
these new models, and the precision study for top pair production in photon-photon collision 
at NLC will be a powerful indirect probe for new physics beyond the standard model. 

This work was supported by the National Natural Science Foundation of China and a 
grant from the State Commission of Science and Technology of China. 
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Appendix A 



The form factors //, f-' are given by 



/I = 



mt 



mi — t 



Vi {Fi - Fo){i, mt, mi) + {Fq - Fi){mj,mt, mi 



=HO,h 



+ ^ r]i\^{Fi + Fo){i,mt,mi) - (Fi + Fo){mt,mt,mi) 

i=A,Z 



i=H+,W+ 



fi = 



{mf - tf 



X ^« ("*t ~ ■ ^^3) fnt, mi) - Fi(m|, mt, mi)^ 

i=HO,h 



+2mlpi ■ p3{Go + Gi){mf,mt, mi) - m] (Fo(t, mt, mi) - FQ(mj,mt, m, 

+ Y m {ml-pi-p3){Fi{i,mt,mi)-Fi{ml,mt,mi)^ 

i=A,Z 



+2mjpi -233(^1 - Go){mf,mt,mi) - mf (^Fo{m^,mt,mi) - Fo{i,mt, 



m. 



+ Y Vi {mt - Pi ■ P3) 0, mi) - Fi(mf , 0, mi)) + 2mtPi • psdimt, 0, 

i=if+,W+ 



mi 



/2 



mt 



t — m 



2 X 'ni\p2 ■ PACui-p2,PA,mi,mt,mt) 



t I i=A,Z 



+ (m^ C21 + pi • P3(Cii - 2Ci2 - 2C23) + 2C24) (Pi, -P3, mi,mt,mt) 
+2 X »?i [P'2 •P4(2Co + Cii)(-]92,]54,mi,mt,mt) 

+ (mf (C21 - 4Co) + pi • P3(2Co + Cii - 2Ci2 - 2C23) + 2C24) (pi, -P3, mi, mt, mt) 

- E [("i?(Co + C21) + pi • P3(Cii - 2C23) + 2C24) (Pl , -p3,mi,mb,mb) 
i=H+,W+ 

-6C24(pi, -P3,mh,mi,mi) +p2 ■ Pa{Co + Cii)(-p2,P4, mj, m;,, m^) 



fV 

/3 — 



t — m 



4 X mfC2i -pi •p3(Ci2 + ^23) {pi,-p3,mi,mt,mt) 



t I i=A,Z 



+4 X »7i ['Ti?(2Cii + C21) -pi •p3(Ci2 + C23) {pi,-p3,mi,mt,mt) 



+2 X { [-?Ti|(Co + Cii + C21) +pi •p3(Ci2 + C23) {pi,-p3,mi,mb,mb) 
i=H+,W+ 



+3 



-mf(Cii + C21) +pi •P3(Ci2 + C23) (pi,-p3,m;„mj,mj)| | 



fV 

/6 — 



i-m^ I i=A,z 



2 ^ rii{2C24:-'rntC2i)i-p2,P4,rrii,'rnt,rnt) 



-2 ^« [2*^24 - m^{ACo + 4Cii + C21 {-P2,P4, m, mt, mt) 



i=HO.h 



fV 

/l2 — 



+ ^ I (Co + Cii + C21) - 2C24 (-P2,P4, mi,mb, nib) 
i=H+,W+ 

+3 2w?t{Cii + C21) - 2p2 • P4(Ci2 + C23) + 2C24 (-P2,P4, mfe, m^, mi)| | 

9 S 'Vi I [2C24 + C21 - 2p2 • P4(Cl2 + C23)l (-f'2,P4, Wi, m*, mt) 

+ [2C24 + m^C2i - 2pi ■ pz{Ci2 + (^23) (pi, -P3,mi, mt, mi)| 

+ XI { 2C24 + mt(C2i - 4Co) - 2p2 •P4(Ci2 + C23) i-P2,P4,mi,mt,mt) 

i=H°,h 

+ [m^ (C21 - 4Co) - 2pi • j03(Ci2 + ^23) + 2C24 {Pi, -P3, mi, nit, "^t)} 



P2 ■ Pi{Ci2 + C23) + -^mtiCo -C21) - C24 



+ 



Pi ■ P3(Ci2 + C23) - -^mjiCo + C21)- C24 



{-P2,P4,mi,mb,mb) 
{Pi,-P3,mi,nib,mb) 



+ ^C24{-p2,P4,mb,mi,mi) + SC24{pi, -P3,mb,mi,nii)^ 



fV 

/l3 — 



mt 



t — mi 



2 'niC2iiPi,-P3,mi,mt,mt) 



[ i=A,Z 

+2 X] Vi{'^Cu+C2i){pi,-p3,mi,mt,mt) 
i=H°,h 

- Y Vi[{Co + Cu + C2i){pi,-p3,mi,mb,mb) 

i=H+,W+ 

+3{Cii + C2i){pi,-P3,mb,mi,mi)] > 



fV 

/it 



mt 



16 



i — mi 



2 X] 'niC2i{-P2,P4,mi,mt,mt) 



+2 XI mC^Cu + C2i)i-p2,P4,mi,mt,mt) 



- X] [(<^n + C2i)(-p2,P4,'rrai,m6,m6) + 3(Cii + C2i)(-i?2,P4,'"^6,?Tii,rni)] 
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where 



-y) + ml{l -y) + mly 



Fn{q,ini,m2) = [ dyy"log 
^0 

Gn{q,ini,m2) = - [ dy — i^—- .,2/1 \ , 2 ' 
7o -q^y{l - y) + mi{l - y) + m^^y 



and C24 = — |A + C24, Co, Cjj are the three-point Feynman integrals, definition and expres- 
sion for which can be found in Ref.[17]. 
The form factor are given by 



M5) 



where 



/r 



6(1) 



f: 



6(1) 



.6(1) 
/3 



6(1) 



,6(1) 
J6 



2mt J2 ?/i(2£'27 - + 2£>3i2) 

i=HO,h 

J2 m [m|(L>o - i^ii + -D12 - D21 - D22 + 2D24 + Dsi - D32 



mt 



=HO,h 



6(1) _ 



-D34 + Dsq) + 2p3 ■ P4(-Dl3 - D12 - £'310 + D26) 

+2(2£>3ii - 3Z>3i2 - 2D27) + 2p2 • P4(^2i - 2£>24 - £'25 + D26 + L'38 - Do) 

+2pi ■ Pa{D22 - -D35) + 2pi • P2(^34 - D^q) 

^ r?i [m|(Di3 + 2D35 + 2D38 - 41)39) - 4p2 • Pi{Di2 - D13 + D23 
i=H°,h 

+D24 - D25 + -D39) - 4p2 • P3(^23 + ^37 + -D310) 
+4ps ■ P4D26 + 4(L»27 - ^312 + 3L>313) 

2 Vi[mU-3Du + 3Di3 + 2D23 + 2D24-2D25-2D26 + D32 + D34 

i=H°,h 

-D35 - D38) + 2p2 • P4{Di2 - Di3 + D24 - D25 - D38 + Dsg) 

+ 2pi ■ P2{D22 + D23 - 2D26 + - Dsio) 

+2p2 ■ P3D37 + 4(D3i2 - -D313) 

2 E Vi[mt{-Dn + D2i + D22+D3i + D3e-Do) 

i=HO,h 

+2pi ■ p4{Du + Di3 + D24 - 2D25 -D35)- 2pi ■ P3{D26 + D310) 

+2pi • P2(^24 + i^34) + 2(2L>27 + i^31l) 

^ J2 Vi [m?(-2r>ii + 2Di3 + D21 - D22 + 2D23 + -Do + 2£>26 - 2£>25 

+2D39 + 2D310) + 2pi ■ P4{-Di2 + Di3 + D25) - 2p2 ■ P4{D23 + ^31o) 
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-2pi-p2D24 - (2L>27 + ^311; 

f^^^^ = Amt Vi{D22-D26 + Ds2-D3e-Ds8 + D3io) 

i=H",h 

f^^^^ = Amt r]i{Du-D^3 + D22-D26-D34 + D35 + D3Q-Duo) 

i=HO,h 

/?^^ = 4 ^ Vi[mt{D24-D26-Ds4 + D3e-Dss + D3io) + D3i3] 

i=H°,h 

/S^^ = 4mt Y Vi{Dn-Di3 + D24-D2e-D3i + D34 + D35-D3io) 

i=HO,h 

fli'^ = -4 E '?i(^27 + ^313) 

i=HO,h 

fi2^ = E r]i[m^t{3Do-2Dn+Di3-D2i-D22-2D23 + 2D25 + 2D2Q + D35 + D38 

i=H°,h 

-D39 + D310) + 2(^27 + 3£>313) + 2pi • P2{-Di2 + Di3 - D23 - D24 + D25 

-D26 - D310) - 2pi ■ P3{D39 + D310) - 2p2 • 

/i3^^ = -2mt E ViiDu + -D22 - D24 + 2D39 + 2D3io){-p2, -pi,Pz, mt,mi,mt,mt) 
i=HO,h 

/li') = 2mt E Vi{D22-D24 + D25-D2e) 

i=H°,h 

f^i^^ = -2mt E Vi{Do + Du-D2i + D24 + 2D39 + 2D3io) 

i=HO,h 

fl^^ = 2mt E m{Do + D12 - Di3 - D21 + D24 + D25 - D2e) 

i=HO,h 

fi?^ = 4 E m{D23-D2e-D3s + D39) 

i=H°,h 

fl8^ = 4 E mi-Dl2 + Di3 + D23-D24 + D25-D2Q + D37-D3lo) 
i=HO,h 

fl9^ = 4 E m{Dl2-Di3 + D23 + D24-D25-D2e + D39-D3lo) 

i=H'>,h 

f20^ = 4 E Vi{D23-D25-D35+D37) 
i=HO,h 

fl^^ = E Vi{'^D312 - D3U) 

i=H + ,W+ 

/2^^^ = E '^i "i-t (-D3I - ^32 - D34 + i:'36) 

i=J?+,W^+ 

+2p2 • P4{Du - D12 + D21 - D24 + £'38) + 2pi • P4(i524 - £»22 - £'35) 

+2pi • P2{D34 - Dse) + 2p3 ■ P4{D2& - D25 - D310) + 2(21)311 - 3i?312) 
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/s^'^ = 2 ^ Vi\mt{D35+D38-2D39-2D3io) + 2p3-p4{D2e-D23) 



10 



+2p2 • P4(-Dl3 - Di2 - -D24 + D25 - -D39 + -D310) 

-2p2 ■ PsDsY - 2pi ■ P2-D310 + 2(^27 - -D312 + 3-D313) 

^ E ["^? (2^24 - 2L'25 - 2L'26 + -D32 + £'34 - -D35 - £'38) 

i=H+,W+ 

+2p2 ■ P4{Di2 - -D13 + -D24 - -D25 - -^38 + D39) + 2p3 ■ P4D23 + 2p2 • P3-D37 
+2pi • P2(^22 - 2L»26 + ^36 - i^31o) + 4(£>312 - £'313) 

^ ^ 77i [m2(D2i + -D22 + D31 + Dae) 

i=H+,W+ 

+2pi ■ P2{D24 + Ds4) - 2pi ■ P2,{D2Q + D310) 

+2pi ■ p4{Di2 - Di3 + D24 - 2L»25 - -D35) + 2(^311 + 2L>27) 

I E ht (^21 - ^22 - 2D25 + 2D2Q + 2D39 + 2D310) 

i=H+,W+ 

+2pi ■ P4(-Dl3 - D12 - D24 + D25) - 2pi ■ P3D310 - (21)27 + D3u) 

mt X! 'ni{D32 - DsQ - Dss + £'310) 

i=H+,W+ 

X! ??j(-D22 - -D24 + -D25 - -D26 - -D34 + -D35 + -D36 - -D310) 
i=H+,W+ 

E rji [mt{D36 - D34 - D38 + -D310) + -D313] 

i=H+,W+ 

X! ?7i(--D21 + £^24 + -D25 - ^26 - ^31 + -D34 + £>35 - -D310) 

i=H+W+ 



i=H+,W+ 

\ E V^[mK-'^Dn + 2Dl3-D2l-D22 + 2D25 + 2D25-D35 + D38 

i=H+,W+ 

-D39 + D310) + 2pi ■ P2{-Di2 + Di3 - D24 + D25 - D26 - -D310) 
-2pi ■ P3(D39 + -D310) - 2p3 • P4D23 - 2p2 ■ P3D37 + 2(D27 + 3-D313) 
^mt E ^i(-^22 + D24 - 2D39 - 2D310) 

i=H+,W+ 

^nit E ^i(-^^22 - D24 + D25 - D2q) 
i=H+,W+ 

Unt E rii{Dii-D^2 + D2i-D24-2D39-2D3w) 

i=H+,W+ 

\mt E Vi{-Dii + D12 - D21 + D24 + £'25 - D2g) 
i=H+,W+ 
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,6(2) 
J17 




V r]i{D23 - D26 - D^s + Ds9) 

i=H+,W+ 




,6(2) 
J 18 




Yl Vi{-Dl2 + Dis + D23 - D24 + D25 - D2Q + r>37 

i=H+,W+ 


- -D310) 


,6(2) 
/l9 


= 


^^(^12 - + D23 + i^24 - -D25 - £>26 + - 

i=H+,W+ 


^310) 


M2) 
J 20 




'ni{D23 - D25 - D35 + D37) 

i=H+,W+ 




rb(3) 
/l 

pb(3) 
J2 




-9mt Y 'ni{-D27 - D311 + D312) 
i=H+,W+ 






,6(3) 

/s 




y 2^ ViK-L'27 + J^312) 
i=H+,W+ 




^6(3) 




-9 ^ m{D3i3 - D312) 

i=H+,W+ 




^6(3) 




9 Y ^^(^27 + -D3I1) 

i=H + ,VV + 




,6(3) 
/6 


= 


-9 Y - ^311) 

i=H+,W+ 




,6(3) 
/7 




-9mt Y ViiDis - D12 - D24. + D25 + D32 - D3Q - 
i=H+,W+ 


D38 + -D310) 


,6(3) 
/8 




-9mt Y ^^(-^13 - D12 - 2L»24 + D25 + r'22 

i=H+,W+ 
-D26 - D34 + + D36 - D310) 




,6(3) 

/g 




-9mt Y ^^(-^13 - - ^21 + ^25 - 2L»34 + ^36 " 

i=H+,W+ 


- D38 + £^310 


,6(3) 
/lO 




-9mt Y ^^(^13 - Dn - 2L>2i + 2L>25 + L»24 

i=H+,W+ 
— D26 — ^31 + — D310) 




,.6(3") 

Jll 

Mi) 

J12 




9 Y '?i^313 

i=H+,W+ 

M3) M3) M3) M3) 
J 13 — J 14 — ^15 — il6 ~ ^' 




,6(3) 

Jii 




-9 ^ rii{D23 - D26 + D^g - Das) 

i=H+,W+ 




,6(3) 
/l8 




-9 Y rii{D23 - D26 + - D310) 
i=H+,W+ 




,6(3) 
/l9 




-9 Y ^^(^23 - D25 + £>39 - -D310) 
i=H+,W+ 




,6(3) 
J 20 




-9 Y ^^(^23 - D25 + £'37 - -D35) 

i=H+,W+ 
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2mt (21^27 - ^311 + 2D312) 

i=A,Z 

mt J2 ["i?(^o - Dri + £>i2 + D21 + D22 + 2r»24 + i^3i - - + Dse) 

i=A,Z 

+2p3 ■ P4{Dis - D12 - D310 + D2%) + 2p2 ■ P4{-Do + D21 - 2D24 - D25 
-D26 + L»38) + 2pi • P4(-2£>13 + £>22 - 2D25 - D35 + 2D12) + 2pi • P2i2D24 

+D34 - Dse) + 2{2D3u - SD^u) 
J2 m [ml{Di3 + 2D35 + 2D38 - 4.DS9) 

i=A,Z 

-Ap2 ■ P4iDl2 - Dl3 + D23 + D24 - D25 + D^g) + 4p3 • ^4-^26 
-4p2 ■ Pz{D2Z + D37 + Dsio) + 4(L»27 - A312 + 3^313 ) 

2 m [mt{-Di2 + - 2D22 - 2D23 + 4L>26 + £'32 + -D34 - D35 - D^s) 

i=A,Z 

+2p2 ■ P4{Dl2 - Di3 + D24 - -D25 - ^38 + ^39) + 2pi • P2(^22 + £>23 
-2L>26 + ^36 - £'310) + 2p2 • P3^37 + 4(L'312 " -D313) 



2 Y 'ni[ml{DQ-Dii+D2i + D22 + D3i + D3Q) 

i=A,Z 

+2pi ■ p4{Di2 + £'13 + -D24 - 2D25 - -D35) + 2pi ■ P2{D24 + D34) 
-2pi ■ P3{D26 + r>3io) + 2(2D27 + -D311) 

2 Y ni[niU-Do-D2i-D22-2D23-2D24 + 4D25 

i=A,Z 

+4L>26 + 2L>39 + 2L>310) + 2pi • P4(-i^l2 + -D13 + D25) - 2p2 ■ P4{D23 + -D310) 
-2pi • P2D24 - {D3U + 2£>27) 

4mt Y Vi{-D22-2D23 + 3D26 + D32-D36-D3s+D3io) 
i=A,Z 

Amt Y '7? (^12 - -D13 + -D22 - D26 - D34 + D35 + Dse - £>3io) 

i=A,Z 

4 Y Vi[mt{-2D23-D24 + 2D25+D26-D34 + D36 

i=A,Z 

-D38 + i:»3io) + -D313] 

4mt Y "niiDii - Di3 + D24 - D26 - D31 + D34 + Das - -D310) 

i=A,Z 

-4 m{D27 + D3i3) 

i=A,Z 
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^ r?i [m2(-3L'i3 - D21 - D22 - 2D23 + + 2D2& + D^^ + L>38 

i=A,Z 

-D39 + D310) - 2pi ■ P3(-D39 + D310) - 2p2 • P3-D37 

+2pi • P2{-Di2 + £>i3 - -D23 - -D24 + D25 - D26 - -D310) 

+2(D27 + 3D313) 

-2mt Vi{-Di2 + 2Di3 + D22-D24 + 2D3Q + 2D3w) 

i=A,Z 

2mt ^ r)i{-2Du + 2Di3 + D22 + 2D23- D24 + D25-3D2e) 

i=A,Z 

-2mt J2 m{Do + Di2- D21 + D24 + 2D3Q + 2D310) 

i=A,Z 

2mt Vi{Do + Di2-Di3-D21+2D23 + D24-D25-D26) 

i=A,Z 

4 Y Vi{D23 - D26 - D38 + D^g) 



i=A,Z 



4 J2 r]ii-Di2 + Di3 + D23- D24 + D25- D26 + - D310) 

i=A,Z 

4 J2 rii{Di2 - Di3 + D23 + D24 - D25 - D26 + Dsg - D31Q) 

i=A,Z 

4 Y ViiD23 - D25 - D35 + D37) 



i=A,Z 



Y m2mtD3i2, = 0, 

i=H+,W+ 

J2 r]im'^D22 - 2D27 - 2p2 ■ P4{D24 + D12 - Di3 + D34) 

i=H+,W+ 

+2p3 ■ P4(£>25 + D35) - P2 ■ P3(2£>26 + ^31o) + 4D311 + m^Dse 
- Yl Vi2D3i3, 

i=H+,W+ 

Y m^i^Sll - D312) + 2p3 ■ P4D35 - 2p2 ■ P4D34 - P2 ■ P3-D3IO 

i=H+,W+ 

+mfD3G + 2p2 ■ P4D36 + 2p2 ■ P3D38 - P3 ■ P4D310 
-m^ D32 + 2p2 • P4(£>22 + D25 - D24 - D26) 

Y ^i2(L>3i2 - £'313), 

i=H+,W+ 

Y r]i2mt{D26 + D310), 

i=H+,W+ 

Y Vi2mt{D3io - D38), 
i=H+,W+ 
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^ r]i2mt{D2e - D22 + -D310 - -Dae), 

i=H+,W+ 

i=H+,W+ 

E ^i2(L>313 - -D311), 

i=H+,W+ 

Yl ViD27, 

i=H+,W+ 

Y mfnt{Di2 + -D24), 
i=H+,W+ 



= 0, 

= Y mm{D2A- D22), 

i=H+,W+ 

= 0, 

rja{D23 - D25 + 1)37 - D»35), 

i=H+,W+ 

Y riaiDsr - D35 - £'39), 

i=H+,W+ 

Y Vi^{D23 + 2D24-2D25-D26 + Di2-Di3 + D34 + D37-D35), 

i=H+,W+ 

Y Vi'^iDu + D37 + D38 - L>35 - -D36 - L>39 + D24 + D26 - D22 -D25), 

i=H+,W+ 

= E r. - Jl V \ {-12mt[m|go - Pz ■ ^4(2^22 - 2C23 + Co) + \]} 

i=HO,h ^ « ' 

(P4, -Pi -]32,mt,mt,mt) + ^ ^{mtCn){-P2,Pi + P2,'mb,mi,mi) 

i=H+,W+ ^* 

/r"" = /8'' = /9^ = /i^o= Y r FT^-12m,[Co + 4(C722 -C23)]} 

(P4, -Pi - P2, mt, mt, mt) 



In the above, Dq, Dij, Dijk are the four-point Feynman integrals [17], and 
Do,Dij,Dijk{-p2,-pi,P3,mt,mi,mt,mt) in f-''^^ and /^'^'^^ 
-Do , Dij , Dijk{-P2 , -Pi ,P3,Q,mi, 0, 0) in /^"^^^ 

Do,Dij,Dijk{-p2,-pi,P3,mi,0,mi,mi) in /f^^^ 
Do,Dij,Dijk{p4,-p2,P3,mb,mb,mi,mi) in Z^'-''' 
and 

S = (pl+p2)^ i=(P3-Pl)^, U={P3-P2f 
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1. 1. 2 ~^ 

P3-P4: = -^S, P1-P4=P2-P3 = ^{mi -U) 

sin^ a cos^ a 

sin^ sm^ p 

2 

r7A = cot^ P, 1]H+ = cot^ /? H 1 tan^ /3, ?7z = = 1, 

Appendix B 

The expressions of Hi{mt,pi ■ p-2^Pi ■ P3^Pi ■ Pa-,P3 ■ Pi,) in the ampUtude squared are given 

by 

Hi = -8mf + Smtpi ■ P2 — SmtPi ■ pz + 8mtP2 ■ P3 

H2 = — 32mf — WmtPi ■ P2 — "i^mtpi ■ ps + 32mtP2 • P3 

H3 = 32m J + 16m^pi • p2 — IGmfpi ■ ps — 8mfp2 ■ Ps 

Hi = 8ml ~ 8m^pi ■ p2 - 16mfpi ■ ps - 16pi ■ P2P1 ■ Ps + 8m^p2 • Pa 

i?5 = 32mfpi • p2 + 16(pi • p2)^ - 16m(P2 • Ps - 16pi • P2P2 • Ps + Smfpi • ps 

Hq = SrutPi ■ p2 - 8ml - 16mfp2 • Ps - 8mfpi ■ ps 

H7 = -8ml + Smfpi • p2 + 4mfpi • p^ - 8mtPi ■ P2P1 ■ Pa + 4m^P2 • Pa 

i?8 = — 8mfpi • p2 + 8mt(pi • ^2)^ + 4mfp2 ■ Pa — ^mfpi ■ p3 

Hg = -8mfpi • p2 + 8ml + 4mfp2 • Pa - 4mfpi • p3 

-H'lo = -8mt{pi ■ P2f + 8mlpi ■ p2 + 8mtPi ■ P2P2 ■ Pa - 4mfpi • ps - 4mfp2 • Pa 
Hii = 8mjpi ■ Pi - 8mjp2 ■ Pa + 8771^^3 • pi - 8pi ■ p4P2 ■ pa 

+8pi • P2Pa • P4 - 8pi ■ pzP2 ■ PA 
H12 = -IGmfpi -^4 - 32mfp2 •P4 + 32m?P3 -^4 - 32pi -pspa •P4 
Hi3 = 16mfpi • p4 + 32mlp2 ■ pa - Sm^Pa • P4 - IGmtPi ■ P3P2 ■ Pa 
Hi4 = —8ml Pi ■ P4 + 8ml p2 ' P4 — 16mtpi • papi ■ p4 + 8mfp3 • p4 

-16mtpi ■ P2P3 • P4 + 16mtpi • ^4^2 • Pa - 16mtpi • P3P2 • P4 
Hi5 = IGmtPi ■ P2P1 • P4 + 32mtpi • P2P2 • P4 - 8mtPi ■ P4P2 ■ Pa + 8mtPi ■ paP2 • P4 
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-Snitpi ■ P2P3 ■ P4 - lGmtP2 ■ P3P2 ■ Pa 
HiQ = -8mfp2 • PA + 8mfpi • p4 - SmtPi ■ P3P2 ■ Pa + 8mtPi ■ P2P2, ■ Pa 

-SmtPi ■ PAP2 ■ Pz - 16mfp3 • p4 
Hv7 = -8m^pi ■ pspi ■ p4 + Amfps ■ p^ + Amfpi ■ p^p2 ■ Pa - 4m^pi • piP2 ■ P3 

+Amlpi ■ P2P3 ■ Pa + Smjpi ■ p4 - 8mjp2 ■ Pa 
His = SnifPi ■ P2P1 ■ pa - Sm^pi • P2P2 ■ Pa - ^mlpi ■ p4P2 ■ pz + 4mf pi • P2P3 ■ Pa 

-4mfpi • P3P2 ■ PA + 8pi ■ P2P\ ■ PAP2 ■ P3 - 8pi ■ P2P1 ■ P3P2 ■ PA + 8(^1 • P2)^P3 ■ PA 

-16pi • P4P1 ■ P2P2 ■ P3 + Sm^pi ■ pspi ■ Pi + 8mlp2 ■ P3P2 ■ Pa - ^rn^ps ■ p^ 
HiQ = 8rr4p2 ■ Pa - Sm^pi • p4 - 4mtPi • pzP2 ■ Pa 

+4mf pi • P2P3 ■ PA - 4mj pi • p4P2 ■ P3 + ^rr4pz ■ Pa 
H20 = Smlpi ■ P2P2 ■ PA - 8mlpi ■ P2P1 ■ Pa - 8mfp2 ■ P3P2 ■ Pa + 4777.^ p3 • p4 

+4mf pi • P4P2 ■ P3 - ^rrifPi ■ p^p2 ■ Pa + 4m|pi • P2P3 ■ Pa 
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Figure Captions 



Fig.l Feynman diagrams contributive to 0{am1 /rriy^) Yukawa corrections to 77 tt: 
(a),(b) tree level diagrams; (c)-(e) self-energy diagrams; (f)-(i) vertex diagrams; (j) in- 
cluding neutral Higgs exchange diagrams; (k) including 77ii'+ii'+(77G+G"'") -coupling 
diagrams; (l)-(n) box diagrams. Here we only plot the one-loop diagrams corresponding 
to tree-level diagram (a). The dashed lines represent H,h,A, H^,G^, for diagrams 
(c),(d),(e),(f),(h),(l), H^,G^ for diagrams (g),(i),(k),(m),(n) and H,h for diagrama(j). 

Fig.2 Plot Aa/ao versus in the standard model. 

Fig. 3 Plot Aa/ao versus for Ma = 600 GeV in the two-Higgs-doublet model {a = j3 = 
0.25). 

Fig.4 Plot Aa/aQ versus Ma for M^, = 600 GeV in the two-Higgs-doublet model {a = (3 = 
0.25). 

Fig.5 Plot Aa/ao versus Ma in the minimal SUSY model for tan/? = 0.25. 
Fig.6 Plot Aa/ao versus Ma in the minimal SUSY model for tan/? = 1. 
Fig.7 Plot Aa/ao versus Ma in the minimal SUSY model for tan/? = 5. 
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